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We present theoretical and experimental considerations pertaining to deeply inelastic heavy- 
flavour production at HERA. The various theoretical uncertainties in the cross section calcula- 
tion are discussed. Cuts are imposed to determine the fraction of charm production accessible 
to the detectors. The production of charm at asymptotic and bottom production are also 
covered. Experimental aspects include current charm production data analysis and prospects 
for future analyses including anticipated high precision and distinguishing photon-gluon fusion 
charm events from excitation from the charm parton density. The feasibility of measuring 
F2^{x,Q^) is investigated. 
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Abstract: We present theoretical and experimental considerations pertaining to 
deeply inelastic heavy-flavour production. The various theoretical uncertainties in 
the cross section calculation are discussed. Cuts are imposed to determine the 
fraction of charm production accessible to the detectors. The production of charm 
at asymptotic and bottom production are also covered. Experimental aspects 
include current charm production data analysis and prospects for future analyses 
including anticipated high precision and distinguishing photon-gluon fusion charm 
events from excitation from the charm parton density. The feasibility of measuring 
F2'^{x,Q^) is investigated. 

1 Introduction 

Heavy-flavour production in deeply inelastic scattering (DIS) at HERA is now emerging as 
a very important means of studying proton structure. The ink is still drying on the first 
experimental reports of charm production from photon- mediated DIS at HERA |1|, 0. The next- 
to-leading order (NLO) calculations have also been published within the last four years. The 
inclusive calculation of the photon-mediated heavy-flavour structure functions F2^{x , Q"^ , m'^) 
[R|, the inclusive sing le differential distribution dF^l/dO [| (O being the transverse momentum 
Pt of the heavy quark and the rapidity y), and the fully differential calculation p are now 
available for a complete NLO analysis of the photon-mediated heavy-flavour structure function. 

In section ^, following a short review of the necessary formulae, we investigate various the- 
oretical issues surrounding DIS heavy-flavour production. The primary sources of theoretical 
uncertainty include the imprecisely determined charm quark mass and renormalization and fac- 
torization scale dependences. Additional impediments to a clean extraction of the gluon density 
from heavy-flavour production include the effects of light-quark (m, s) initiated heavy-flavour 
production and the influence of the longitudinal heavy-flavour structure function Q^, m?) 

upon the cross section results. 

We investigate the effect of realistic cuts in pt and pseudorapidity rj on the cross section 
and determine acceptance probabilities as a function of x and Q^. Charm production in the 



limit ^ ml and the transition of charm production from boson-gluon fusion at low 
to excitation from the charm density as becomes much larger than is then discussed, 
followed by a cursory view of bottom production. 

In section ^, we discuss the analysis of the 1994 HERA data. Identification methods of 
the produced D^,D*^ are outlined. The measured cross section is compared with various 
theoretical predictions. A determination of the source of the charm production is performed, 
revealing the primary production mechanism at presently measured values is photon-gluon 
fusion rather than stemming from the charm parton density. The charm structure function F2'^ 
is extracted and the ratio F2'^/F2 is determined. 

Future experimental prospects include the installation of silicon vertex detectors, enabling 
greater charm and bottom hadron detection efficiency. The anticipated luminosity of 500 pb~^ 
will allow detailed studies of charm production dynamics. The transition from boson-gluon 
fusion of charm to excitation from the charm quark sea should become apparent as the accessible 
Q2 grows. The predicted bottom quark production cross section will enable studies of F2^/ Fl'^ 
as a function of x and with reasonable precision. 



2 Theoretical Aspects 

2.1 Background 

The reaction under study is 

e-(/) + Pip) ^ e-(/') + g(pi)(g(pi)) + X ^ 



(1) 



where P{p) is a proton with momentum p, Q{pi){Q{pi)) is a heavy (anti)-quark with momentum 
Pi (Pi = ''^^) ci-nd X is any hadronic state allowed by quantum number conservation. The cross 
section may be expressed as 



dxdQ^ xQ^ 



where 



q = l-l', Q' = -q 



1 + (1 - Q^ m') - y'Fl\x, Q\ m') 



X 



y 



p-q 



2p ■ q ' " p ■ I 

The inclusive structure functions F^^ were calculated to next-to-leading order (NLO) in Ref. 
The results can be written as 



(2) 



(3) 



Fk{x,Q\m^) 



+ 



A'K'^m? Jx z 



,x 



■{ei/.(f,/i^)(cS + ctJln^ 



'=9,9 



(4) 



where k = 2,L and the upper boundary on the integration is given by 2;max = Q^/ {Q^ + 4m^). 



The functions fi{x,jj'^ 



g, q, q) denote the parton densities in the proton and /i stands 



for the mass factorization scale which has been set equal to the renormalization scale. The 

cg(C,0(« = 9,Q,r.i = 0,1), cg(C,0(^ = 9,q,q-J_= i), and <kC,0(^ = Q^rJ = i) 

are coefficient functions and are represented in the MS scheme. They depend on the scaling 
variables ( and ^ defined by 

c = A-i ^ = —2- (5) 

where s is the square of the cm. energy of the virtual photon-parton subprocess Q^{1 — z)/z. 
In Eq. we distinguish between the coefficient functions with respect to their origin. The 
coefficient functions indicated by cl^\{C,0^'^k\iC^O originate from the partonic subprocesses 
where the virtual photon is coupled to the heavy quark, whereas dl^\{C-,0 comes from the 
subprocess where the virtual photon interacts with the light quark. The former are multiplied 
by the charge squared of the heavy quark efj, and the latter by the charge squared of the light 
quark e| respectively (both in units of e). Terms proportional to enei integrate to zero for the 
inclusive structure functions. Furthermore we have isolated the factorization scale dependent 
logarithm ln(/i^/m^). 

The scale dependence and the poorly known charm quark mass are the largest contributors 
to the theoretical uncertainty. The effects of F^"^ and the light-quark initiated contributions 
are discovered also to be important in the analysis. To aid the experimental analysis, the fully 
differential program p is used to apply a series of cuts to determine the percentage of events 
the detectors are likely to see in bins of x and Q^. With the planned inclusion of silicon vertex 
detectors, the ability to see bottom events increases dramatically, motivating the presentation 
of results for the cross section and F2^{x, Q^, ml). HERA is in a unique position to evaluate the 
transition of charm production from photon-gluon fusion to excitation from the charm parton 



density. For other phenomenological investigations, see H, 0, |T0|. 



2.2 Code Update 

For this study we use an updated version of the code based upon The original code was 
based upon fitting the coefficient functions described in eq. using a two-dimensional tabular 
array of points in ( and ^. The coefficient functions were generated via a linear interpolation 
between the calculated points. The linear interpolation was insufficiently accurate and required 
a more sophisticated interpolation procedure. 

The present interpolation procedure is based upon a Lagrange three-point interpolation 
formula, see eq. 25.2.11 [|ll[]. The results have been thoroughly compared with the original 



code in and [^] and excellent agreement has been established.^ 

As a demonstration of the code, we present results for the Born C2°''(C,C = 1, 10) in Fig. |l] 
g(C,e= 1,10) and c2( 



and 4;j(C, e = 1, 10) and c^^J (C, ^ = 1, 10) in Fig. |. 



2.3 Scale and Parton Density Related Issues 

The most important numerical sources of theoretical uncertainty in DIS heavy-flavour produc- 
tion are the factorization/renormalization scale dependence and the poorly known charm quark 



^The code is available at |http : / / www. ifh .de/thcory/ publist . html 
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Figure 2: The NLO coefficient functions C2^g{C,,^), for ^ = 1, 10 (upper curves at large C,) and 
C2^g(C)0 fof C = 1) 10 (lower curves at large C,). 



mass. Varying /i in Eq. (Q) indicates the stability of the NLO result against scale changes. To 
get a concrete idea of the effect of scale variations, we first construct a "data" set. The number 
of DIS charm events is calculated for an integrated luminosity of 500 pb~^, using CTEQ3M 
p!2| parton densities with A4 = 239 MeV, rric = 1.5 GeV, and /x^ = + 4m^. Unless oth- 
erwise mentioned, CTEQ3M is used for all results. The results in Fig. |^. are produced for 
1.8 < < 1000 GeV^ and 10~^ < x < 1 using four bins per decade for both and x. The 




Figure 3: Projected number of DIS charm events for an integrated luminosity of 500 ph ^ binned 
in X and with no cuts applied. 

"data" are peaked strongly at small x and Q^. As grows, the events become more evenly 
distributed in x. For small Q^, the number of events falls off as x — > 1. At larger Q^, N 
rises, peaks at intermediate x, then drops. 

Using the "data" set as a point of reference, we can investigate where in the kinematic 
region the effects of scale variation are most strongly felt. Keeping every other parameter fixed, 
we determine the number of events for /i^ = 4m^ and /i^ = 4((5^ + 4m^). We investigate the 
quantity 

Ar(/i2 = Ami) - N{fi^ = 4(Q2 + 4ml)) 



^,{x,Q')^ 



(6) 



2A^(/i2 = Q2 _^ Ami) 

The results are displayed in Fig.. ^ Apart from the high-x region, which contains very few 
charm events, the scale dependence varies relatively little with Q^. As x 0, we find the scale 
dependence disappearing. This behaviour bodes very well for a low-x extraction of the gluon 
density from charm events. 

The charm quark mass uncertainty presents a stickier problem. A precise measurement of 
the charm mass is yet to be made. To develop a feeling for how much of an effect the uncertainty 
has, we calculated 

, . ^2, _ Njm, = 1.3GeV) - Njm, = 1.7GeV) 

2N{m, = 1.5GeV) ' 



Very large effects are naturally found near threshold, but as Q'^ increases and x decreases, A^^ 
approaches a value on the order of 0.1. Varying the charm mass from 1.3 to 1.7 GeV is a 



conservative estimate; the error induced by the uncertainty viewed in Fig. ^. can be viewed as 
an upper bound. 

A clear indication of the abihty to extract the gluon density from charm production is 
whether one can distinguish the gluon densities from different available parton densities. We 
compare the cross section generated with CTEQ2MF[|l2| (with a flat gluon density as a; — > 0) 
with GRV94HO[0 (with a steep gluon density as x —>■ 0). We define 

A'CTEQ2MF — -^GRV94HO ^g-^ 
2 A/'cTEQSM 

and show the results in Fig. ^. Away from large x, Agiue(x,(5^) is flat as a function of Q^. 
In the intermediate region in x, very little distinguishing power is observed. Beginning near 
X = 10~^, a definite difference is seen. The analysis must seemingly extend to x < 5 ■ 10~^ to 
distinguish cleanly the gluon densities mentioned. 

To summarize: While the scale dependence is well under control, the dominant source of 
uncertainty is clearly the charm quark mass, which has a strong influence on the cross section 
at low X and low Q^. This region is exactly the region sensitive to the gluon density. This 
strong influence of rric poses problems for a clean extraction of the gluon density at small x 
from inclusive measurements. 



Agiue(a;,Q^) = 



2.4 Smaller Contributions to the Cross Section: 
Quark Initiated Results 



pec 



and Light- 



A clean extraction of the gluon density may be hindered by contributions of Fl'^ and light-quark 
initiated contributions to the cross section. We investigate the fractional Fl contribution to 
the cross section by plotting 



C 



Fl 



(9) 




Figure 4: A^(x,Q^) (see Eq. (^) binned in x and . 



Figure 5: Am^(x, Q^) (see Eq. ^) binned in x and . 




Figure 6: Agiuc(x,(5^) (see Eq. ^) binned in x and . 

where represents the contribution to the cross section in Eq. (Q) with set to zero. The 
results are shown in Fig. |^. We find a sizeable contribution to the cross section at high y. 
This overlaps with the low and low x region previously determined to be deemed the most 
suitable for a gluon density extraction. To do so, however, one must take into consideration 

We investigate the light-quark initiated contribution by determining 

(10) 

and displaying the results in Fig. ^. The results are very nearly constant, amounting to a 5 - 8 
% contribution except at large x, where the charm contribution does not appreciably contribute 
numerically. 
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Figure 7: Ci?^(x, Q^) (see Eq. ^) binned in x and . 




Figure 8: Ciq (see Eq. ^Tdj)) binned in x and Q^. 

Summarizing: the contributions to the cross section from F^'^ are noticeable in the large y 
region, which overlaps with the small x and small region, hindering the extraction of the 
gluon using only F2'^. The light-quark initiated cross section contributes on the order of 5 % to 
the total cross section, and therefore cannot be totally neglected. 

2.5 Cuts vs. No Cuts 

With the fully differential code developed in a series of cuts can be applied determining a 
more realistic expectation of charm events able to be detected. 

The cuts imposed upon the data are 



(11) 







a (nb) 


(^cuts (nb) 


Efficiency (%) 


O.lOOOE-03 


0.1778E-03 


0.439E+01 


0.106E+01 


24.1 


0.1778E-03 


0.3162E-03 


0.414E+01 


0.113E+01 


27.3 


0.3162E-03 


0.5623E-03 


0.360E+01 


0.961E+00 


26.7 


U.50/3E-U3 


U.1UUUE-U2 


U.29oE+Ul 


U.7U4E+(JU 


23.8 


O.lOOOE-02 


0.1778E-02 


0.232E+01 


0.396E+00 


17.1 


0.1778E-02 


0.3162E-02 


0.174E+01 


0.171E+00 


9.8 


0.3162E-02 


0.5623E-02 


0.125E+01 


0.381E-01 


3.0 


0.5623E-02 


O.lOOOE-01 


0.853E+00 


0.189E-02 


0.22 



Table 1: 
GeV. 



Cross sections with and without the cuts mentioned in Eq. ([IJj for 2 < < 10 





•^max 


a (nb) 


Gcuts (nb) 


Efficiency (%) 


O.lOOOE-03 


0.1778E-03 


0.334E+00 


0.584E-01 


17.5 


0.1778E-03 


0.3162E-03 


0.116E+01 


0.251E+00 


21.6 


0.3162E-03 


0.5623E-03 


0.173E+01 


0.484E+00 


28.0 


0.5623E-03 


O.lOOOE-02 


0.198E+01 


0.681E+00 


34.4 


O.lOOOE-02 


0.1778E-02 


0.191E+01 


0.768E+00 


40.2 


0.1778E-02 


0.3162E-02 


0.162E+01 


0.701E+00 


43.3 


0.3162E-02 


0.5623E-02 


0.127E+01 


0.504E+00 


39.7 


0.5623E-02 


O.lOOOE-01 


0.945E+00 


0.280E+00 


29.6 


O.lOOOE-01 


0.1778E-01 


0.660E+00 


0.112E+00 


17.0 


0.1778E-01 


0.3162E-01 


0.430E+00 


0.254E-01 


5.9 


0.3162E-01 


0.5623E-01 


0.256E+00 


0.136E-02 


5.3 


0.5623E-01 


O.lOOOE+00 


0.131E+00 


0.468E-05 


0.0036 



Table 2: 
GeV. 



Cross sections with and without the cuts mentioned in Eq. (fllj for 10 < < 100 



rjc being the pseudorapidity of the detected charm quark. Looking first at the low range 
from 2 to 10 GeV^, we present the cross section with and without the aforementioned cuts in 
the X bins outlined in Tab. [1|. These results were generated using GRV94HO, rric = 1.5 GeV, 
3 light flavours, and A3 = 248 MeV. In the low x range, we observe an efficiency of 20 - 25 %, 
diminishing near threshold, the low efficiency mostly a result of the pt cut. 

Table 1^ displays results in the medium range, from 10 to 100 GeV^. We find the efficiency 
has risen considerably to the 30-40 % range where the cross section is peaked in intermediate 
X, with on the order of 20 % for the small x region. 

In the large range, from 100 to 1000 GeV^, the results are found in Tab. ^ The efficiency 
continues to rise with Q^, reaching 50-60 % in the high to intermediate x region and 30-40 % 
in the small x range. 

The number of events landing within the imposed cuts mentioned in Eq. 



11) is not over- 



whelmingly large, but a large enough sample should be accessible to gather significant statistics. 







a (nb) 


CTcuts (nb) 


Efficiency (%) 


O.lOOOE-02 


0.1778E-02 


0.401E-01 


0.120E-01 


29.9 


0.1778E-02 


0.3162E-02 


0.124E+00 


0.486E-01 


39.2 


0.3162E-02 


0.5623E-02 


0.159E+00 


0.885E-01 


55.7 


0.5623E-02 


O.lOOOE-01 


0.157E+00 


0.958E-01 


61.0 


O.lOOOE-01 


0.1778E-01 


0.128E+00 


0.757E-01 


59.1 


0.1778E-01 


0.3162E-01 


0.906E-01 


0.458E-01 


50.6 


U.31d2E-U1 


U.5o23E-Ul 


U.572E-(J1 


U.173E-U1 


30.2 


0.5623E-01 


O.lOOOE+00 


0.315E-01 


0.285E-02 


9.1 


O.lOOOE+00 


0.1778E+00 


0.140E-01 


0.167E-03 


1.2 


0.1778E+00 


0.3162E+00 


0.407E-02 


0.400E-05 


0.098 


0.3162E+00 


0.5623E+00 


0.466E-03 


0.413E-08 


0.0009 



Table 3: Cross sections with and without the cuts mentioned in Eq. ([71[) for 100 < < 1000 
GeV. 

2.6 Charm Production at Asymptotic 

Near the threshold for charm production the deep inelastic structure functions Fj, i = 2,L, 
which include the contributions of the light partons u, d, s, and g and the charm quark with 
mass nic are given by 

F.(g^ ml 3) = I [E(/.^ 3) ® {Cl{%, s) + ^, 3) }] 



9 



■)2 ^2 



+G(^^3)®{C|;,(^.3) + Lf„(|J.^,3)} 

+A(^^3)»{c,■::fg.3)+L■^f(^.^;|,3)} 
+i[E(^.^3)®/^-(g.^|.3; 
+G(,.^3)®i/.?i-'2' 



(12) 



where we have the following definitions. The variable denotes the virtual mass squared 
of the photon exchanged between the electron and the proton.. The momenta of the photon 
and the proton are given by q and p respectively and the Bjorken scaling variable is defined 
by X = Q^/{2p ■ q) (g^ = —Q^ < 0). The structure functions Fj depend on x, on and on 
the charm quark mass mc- The convolution symbol ® appears on the right-hand-side which is 
defined by 



{f®g){x) 



1 rl 

dzi dz2 6{x - Z1Z2) f{zi) g{z2) . 
Jo 



(13) 



The gluon density is denoted by G(/i^, 3) where the number of light flavours is taken to be 
three. Also the singlet and nonsinglet combinations of the light quark densities are defined for 
n/ = 3 as follows 

S(/i^ 3) = M(/i^ 3) + M(/i^ 3) + c/(/i^ 3) + J(/x^ 3) 

+s(/i^3) + s(/i^3), (14) 



2 



9 



A(/i^ 3) = - u{fi\ 3) + u{i2\ 3) - - d{n\ 3) + d{iJ,\ 3) 



9 



+s(/i^3)+s(/i^3)J. (15) 

The same singlet and non-singlet classification can also be made for the light parton coefficient 
functions Ci^k {i = '2,L,k = q,g) and the heavy quark coefficient functions Lj^, Hi^, Hi ^ 
depending on rric- The coefficient functions Lj ^, Hi k reflect the corresponding production 
processes. The functions Hi^k contain the processes where the virtual photon couples to the 
heavy quark (c./. Figs. 2, 3, 4, 5. a and 5.b in 0]), whereas Li^k describe the reactions where 
the virtual photon couples to the light quark (c./. Figs. 5.c and 5.d in 0). Hence Li^k and Hi^k 
are multiplied by ef (where I]f=i = 2/9) and by el (where Cc = 2/3) respectively. These 
charge factors are exhibited by Eq. (0). Furthermore both the coefficient functions and the 
parton densities depend on the mass factorization scale n , which, for convenience will be set 
equal to the renormalization scale. The latter also shows up in the running coupling constant 
asifi'^, 3) where the number of light flavours is three. 

Equation (^2|) gives an adequate prescription as long as the cm. energy is not too far above 
the charm threshold, which implies that is not too large compared to m^. However, when 
we enter the asymptotic region ^ , the heavy quark coefficient functions behave like 
ln*(m^//i^) ln-'(Q^/m^) so that the higher order corrections can become large. At sufficiently 
large the charm quark should be treated in the same way as the light partons were at smaller 
Q^. The logarithmic behaviour of the coefficient functions is due to the collinear singularities 
which are regulated by rric. Therefore when ^ m^, the charm quark behaves like a massless 
quark similar to the behaviour of the normal light quarks {u, d, and s) over the whole range. 
Following the same procedure as has been used for the light partons, the mass singular {rric 
-dependent) terms have to be factorized out of the heavy quark coefficient functions using the 
method of mass factorization. This leads to a redefinition of the parton densities in Eqs. (|I^ , 
and the heavy quark coefficient functions turn into the light parton analogues, wherein the 
number of light flavours is enhanced by one. The above procedure is called the variable flavour 
number scheme (VFNS) which is outlined in leading order in |T^ . 



Since all coefficient functions are now available up to order a^, this analysis can be extended 
to NLO to give a better description for the structure functions Fi{x, Q^) at large Q^. A preprint 
is in preparation [TH]. 



2.7 Bottom Quark Production 

The anticipated amount of bottom production is greatly reduced due to the reduction of the 
charge factor by four as well as a significantly reduced phase space. Given a luminosity of 500 
pb~^, however, thousands of events are expected per bin of x and as shown in Fig. ^. We take 
nib = 4.75 GeV, the factorization scale fi"^ = Q"^ + m^, as well as the CTEQ3M distributions. 
A discussion of the implications of this number of b events can be found in section p.2.2| . 



3 Experimental Aspects 




Figure 9: Projected number of DIS bottom events for an integrated luminosity of 500 pb ^ 
binned in x and with no cuts applied. 

3.1 Current Experimental Situation 

Results on charm production in deep inelastic ep scattering are available from both experiments 
HI and ZEUS based on an integrated luminosity of roughly 3 ph~^ collected with each exper- 
iment at HERA in 1994. The HI collaboration has performed the tagging of heavy quark 
events by reconstructing D'^(1864)0 and D*+(2010) mesons, while the ZEUS collaboration 
has given preliminary results for the inclusive D*+(2010) analysis. 

The is identified via its decay mode —>■ K~7c'^ and the D*^ through the decay chain 
D*^ D^Ti^i^^ JC^TT+vr^Q^. For the latter use is made of the tight kinematic constraint 
for the decay of D*^ — > -D^vrJo^ A better resolution is expected for the mass difference 
Am = m{D^iTfi^^) — m{D^) than for the D*^ mass itself. 

Figure [1^ shows the mx-K distribution obtained in the inclusive analysis of HI and the 
Am distribution as observed in the inclusive D*^ of ZEUS. Evidently the number of observed 
events containing heavy quarks is small. Only of the order of 100 to 200 charm mesons are 
identified in any of the different analyses. Combining the and D*^ analysis of HI leads to 
a charm production cross section of 

a{ep eccX) = 17 A ± 1.6 ± 1.7 ± 1.4 nb (16) 

in the kinematic range 10 GeV^ < < 100 GeV^ and 0.01 < y < 0.7. The errors refer to the 
statistical, the experimental systematic and the model dependent error, respectively. The model 
dependent error accounts for the uncertainty in the determination of acceptances due to the 
choice of the parton density in the proton, the mass of the charm quark, and the fragmentation 
function. This cross section is somewhat larger than predicted by the NLO calculations 0, 
The prediction from the gluon density in the proton extracted from the NLO QCD fit of HI [|17 
comes closest to the charm data. 



^Charge conjugate states are always included. 
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Figure 10: (a) The K^tt^ mass distribution observed in the inclusive analysis of DIS events 
from HI, and (b) the Am distribution obtained in the inclusive D*^ analysis of DIS events 
from ZEUS. 
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Figure 11: Normalized xd distribution in deep inelastic ep scattering at (W) ~ 125 GeV for 
\r]D\ < 1-5. The open/closed points represent the /D* data of HI. The shaded histogram 
shows the BGF expectation according to AROMA. The dashed histogram shows the charm sea 
expectation, obtained by selecting QPM events from LEPTO. The full line gives the result of 



(-) 



the QCD evolution of the fit to the charged current v N data 



Information on the charm production mechanism in neutral current DIS at HERA is ob- 
tained from the distribution of xc = 2\P^o\/W , where P^o denotes the momentum of the D'^ 
in the Yp system. Figure |ll| shows the distributions 1/cr da/dxn, which are a convolution 
of the charm production spectrum with the fragmentation function compared to the LO BGF 
expectation of the AROMA generator. The BGF model agrees very well with the shape of 
the data. The figure also includes the expectations for charm mesons originating from quarks 
in the proton either by using LEPTO 6.1 generator pOl, from which only charm sea quark 



events are selected, or by extrapolating the results from charm production in charged current 
^1/ N scattering [21] to HERA energies. Large differences in the shape of the distributions are 
observed for the data and these QPM expectations. From these differences, it is concluded 
that more than 95% of the charm production in neutral current deep inelastic ep scattering is 
due to boson-gluon fusion. This observation seems to be in contradiction to recent inclusive 
calculations of charm production in DIS |T^, from which was concluded that for the kinematic 
range accessible in the current analyses at HERA charm quarks may already be considered as 
partons in the proton [ITU . 



The measurements of the charm contribution -F2^(x, Q^) to the proton structure function 



• H1 D' o H1 D° □ EMC 




Figure 12: F^^^ as derived from the inclusive D*^ (full dots) and analysis (open circles) in 
comparison with the NLO calculations based on GRV-HO (full line), MRSH (dashed line), and 
MRSDff (dash-dotted line) parton distributions using a charm quark mass of m^. = 1.5 GeV for 
(Q^) = 12, 25 and 45 GeV^. The inner (outer) error bars refer to the statistical (total) errors. 
The shaded band represents the prediction from the HI NLO fit to the F2 measurements. The 
EMC data are also shown (open boxes). 
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with decreasing x. 
the MRSH |2|], 



and the MRSDO' ||25| parameterizations of the gluon density in the proton for a charm quark 
mass rric =1.5 GeV. The data are also compared to the prediction from the HI QCD fit to 
the F2 measurements using a charm quark mass of rric = 1.5 GeV. The error band shown for 
this fit includes the propagation of the statistical and the uncorrelated systematic errors on the 
total F2 data through the fitting procedure. This prediction lies systematically above all other 
calculations, independently of x and but agrees better with the ^2^*^ measurement. Within 
the errors the ratio of (^F2^/F2^ is found to be independent of x and Q^. Averaged over the 

kinematic range a ratio (^F2'^ / F2^ = 0.237 ± 0.021 ± 0.041 is obtained, which is one order of 
magnitude larger than at larger x measured by EMC. 



3.2 Future 

3.2.1 Charm Production 

In this section some estimates about the precision that may be expected at HERA for an 
integrated luminosity of 500 pb~^ are given. As an example the capabilities of the HI detector 



26| are considered. For heavy flavor physics an important new feature of the apparatus is the 
double layer sihcon vertex detector (CST) ^\ which will allow use of the apparent proper time 
of charm hadrons in selecting heavy flavor events.. Other decay channels than —>■ K^-k^ 
and D*^ D^ufi^^ — > K'n'^nfi^^ are already under investigation or will become accessible 
due to the CST. 

Table ^ contains a list of decay modes which should be feasible for charm tagging in HI in 
future by also including the CST. A total charm selection efficiency of 1% may be obtained. 
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0.06 
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0.091 


0.07 


0.0016 
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0.07 


0.07 


0.0012 


Sum 


CST 


0.25 


0.170 




0.0031 


Sum D 










0.0097 



Table 4: Compilation of various decay channels of charm mesons accessible in HI. The 
channels marked by "CST" will be accessible only by using the silicon vertex detector. BR{D — *■ 
FS) denotes the product of all branching ratios involved in the decay into the final state "FS". 
For the determination of the total charm tagging efficiency the correlations in the and D*^ 
analyses are taken into account. 



Compared to the present analysis this corresponds to an increase in the total selection efficiency 
of a factor 4 to 5. Although the new channels opened by the use of the CST corresponds only to 
50% of this gain, the signal to background ratio will improve for all decay modes considerably. 
Because of this the effective number of events per luminosity iVg// = {N signal I CFstatY will increase 
by a factor of 6 to 8. Due to the cut in the impact parameter of the mesons, which is 
likely to be different in the and D*^ analyses, the events selected in the D*^ analyses 
will only partially overlap with the events selected in the corresponding analyses. Taking 
an integrated luminosity of 500 pb~^ about 160,000 tagged charm events are expected for 
the kinematic range of the published HI analysis. This has to be compared to the currently 
analyzed 250 events. Taking into account the increase in the effective number of events a gain 
of a factor of 1000 is expected for the statistical significance. In total we expect to observe 
roughly 6,000 double tag charm events in the range 1.7 GeV^ < g2 < 560 GeV^ which would 
allow a study of the charm production dynamics in detail (see Ref. [0)- 

Figure |13| shows the result of the hypothetical measurement of F^^^ in the range 1.7 GeV^ < 
< 560 GeV for a luminosity of 500 pb~^ based on the gluon density determination from 
the NLO HI fit to the inclusive F2 data combining the statistics of the D meson decay modes 
summarized in Table § by assuming rric = 1.5 GeV and /i^ = + 4m^. The contribution due 
to Fl and the light flavors are not included. Statistical and full errors, obtained by adding 
the statistical and experimental systematic error in quadrature are shown. The change in the 
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Figure 13: Expected F^^^ for a luminosity of 500 pb^^ . The points show the prediction from 
the gluon density determination by the NLO HI fit to the inclusive F2. The inner (outer) 
error represents the statistical and the total experimental error. The full (dashed) line gives 
the expectation from the NLO calculations based on GRV-HO and MRSH parton distributions 
using a charm quark mass of = 1-5 GeV. 



acceptance due to the value of rric and /i is not included in the error. For < 100 GeV^ the 
precision of this measurement will be limited by the experimental systematic uncertainty. 

A detailed analysis of the experimental systematic uncertainties in the determination of F2'^ 
with the current experimental statistics can be found in Ref. |]l[ for the HI analysis. At present 
the total systematic error is approxiamtely 20 %. It is dominated by the uncertainty in the 
assumptions made to extract the signal from the observed mass distributions. This is clearly 
an effect of the limited statistics available in the current analyses. Because of the increase in 
statistical significance of the data with a luminosity of 500 pb~^ a much better understanding of 
the measured mass distributions as well as of the detector is expected. Ultimately a systematic 
error of 10% is achievable, which will then be to equal parts due to detector and analysis related 
errors (7%) and to the knowledge of the fragmentation probability P{c ^ D) of a. charm quark 
into a specific charm meson and their branching ratios (7%). For the latter the experimental 
situation is not expected to improve in the near future. 

Figure ^ also shows the NLO predictions using the GRV-HO (1992) and the MRSH param- 
eterization of the gluon density in the proton using = 1.5 GeV and /i^ = + 4m^. For a 
given value of rric the data will still allow a sensitive indirect determination of the gluon density 
even with the relatively large experimental systematic uncertainties. Due to the high statis- 
tics it will be possible to measure the charm production cross section up to > 1000 GeV^. 
Assuming that charm tagging is performed for all decay modes listed in Tab. ^ the sensitivity 
limit for the inclusive measurements will already be reached at 50 pb~^ for < 100 GeV^, at 
which time the experimental systematic will dominate in this kinematic range. 

The discussion in the theory section has shown, that the measurement of Fg'^ is sensitive 
to the gluon density to some extent at large y. Unfortunately the inclusive measurement in 
this range is also very sensitive on the charm quark mass. Therefore it is questionable whether 
the measurement of alone will allow an extraction of the gluon density at small x. In the 
experimental analysis of Ref. it was observed that the change in F2'^ at small due to 
rric is compensated by the change in the acceptance such that the measurement should still 
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Figure 14: Relative change Axd in the shape of the distribution 1/a da/dx£) according to the 
LO Monte Carlo Simulation of the AROMA program for 6 GeF^ < g2 < iqo GeV^ and the 



cuts of Eq. (\n\). The full line gives the change in the shape by changing from 1.3 GeV 
to 1.7 GeV. The dashed line shows the influence of using the GRV-HO (1992) instead of the 
MRSH parameterization of the gluon density in the proton. 



allow the extraction of the gluon density from the inclusive measurement. The behavior of 
the acceptance at < 10 GeV^ has not been studied yet. If it happens that F^^^ does not 
provide an reliable indirect extraction of the gluon density because of the uncertainty in rric, 
exclusive distributions of the identified charm hadrons have to be studied. As an example the 
result of the Monte Carlo study of the influence of rric and the gluon density on the distribution 
1/a da/dxD is shown in Fig. H for 6 GeV^ < < 100 GeV^ and the cuts of Eq. (0). Here 



/S. xd IS defined as 

l/aidai/dxD{xD) 
^^D = 77 — 1 — — 7 — r 17 

where i,j denote the different values of the parameters, i.e. rric = 1-3, 1.7 GeV and PDF=MRSH, 
GRV-HO(1992), respectively. The mass of the charm quark affects strongly the shape of the 
xd distribution (full line) while there is only very little effect due to the choice of the par- 
ton density function (dashed line). A study of this distribution, for instance, should therefore 
disentangle the effect of the charm quark mass and the gluon density on Fg*^. 

In section 2.5, the question when charm quarks should be treated as partons of the proton 
was discussed. The investigation of the x and dependence of the exclusive measurement of 
1/a da/dxD (see Fig. 0) will allow a study of the origin of charm production as a function of 
the kinematic variables. At sufficiently large it is expected that the charm quark behaves 
like a parton in the proton. This will result in a change of the xo distribution from the boson 
gluon fusion dominated regime presently observed in the available data at (Q^) ~ 25 GeV^ to 
the sea quark dominated regime at large Q^. If it becomes possible to control the /z and rric 
dependence of charm production at large Q^, the measurement of F|'^(a;, Q^, xd) can determine 
g{x, Q"^) and c{x, Q'^) simultaneously. This study would certainly require a luminosity of 500 
ph~^ to produce enough data at large Q^. 



3.2.2 Bottom Production 

In section 2.6, predictions for the number of bottom events as a function of x and were 
given for a luminosity of 500 pb~^. In the following the possibilities to measure F2^{x,Q'^) at 
HERA will be discussed using again the HI detector as an example. 



Compared to charm quark events the major experimental difference in bottom quark pro- 
duction is the relatively long lifetime of the B mesons. With use of the CST, bottom events 
are selected by applying a cut in the impact parameter of tracks not fitting to the primary 
event vertex. The combinatorial background is negligible for this analysis, while the charm 
production is a significant background source, because of the much larger cross section. Two 
different possibilities will be discussed here, namely 

1. the exclusive analysis of reconstructed D mesons. This method benefits from 

(a) the branching of B mesons into charm meson being very large [^, i.e. BR{B — >■ 
D^X) = 0.242 ± 0.033, BR{B D^/D°X) = 0.58 ± 0.05 and BR{B D*^)X = 
0.0231 ± 0.033, and 

(b) the relatively long visible lifetime observed in these decay chains. Taking into account 
the decay length of both the B mesons and the subsequent D mesons, values of 
CT ~ 590 fim and cr ~ 780 /im are obtained for the decay chains B —>■ D^/D X 
and B D^X, respectively. This has to be compared with the decay length the 
charm mesons of cr = 124 /im for the / D and of cr = 317 /im for the D^. 
Requiring a selection efficiency for the cut in the impact parameter of 50% would 
reduce the contamination of charm events by a factor of 12.5(2.8) in case of the 
D°{D^) analysis. 

All decay modes summarized in Tab. ^ will also be accessible in this case. Due to the 
larger cut values in the impact parameter for the different decay modes the combinatorial 
background is smaller. This will enable us to use softer cuts. Compared to the measure- 
ment of the charm production cross section in total an increase in efficiency by a factor 
3 may be expected. A systematic error of 16% may be achieved. It will be dominated by 
the uncertainties in the branching fraction BR{B DX). 

2. the analysis of the impact parameter distribution of at least 3 to 4 tracks fitted to a 
common secondary vertex. Compared to the exclusive analysis this method leads to a 
higher selection efficiency, but will be limited by larger systematic uncertainties in the 
subtraction of the charm quark induced background. 

Figure |l^ shows the ratio of ¥2^^ / F2'^ expected for the exclusive D meson analysis from 
bottom quark production using the the calculations of sections 2.1 and 2.6. The charm pro- 
duction background is subtracted statistically. The errors refer to the statistical and the total 
experimental error. For a luminosity of 500 pb~^ the error on the measurement of this ratio 
will still be statistics limited in most of the x and plane. For a given or x, this ratio is 
expected to rise with decreasing x or increasing Q^, respectively. Integrated over the kinematic 
range a mean value of F|7Ff ^ 0.02 is predicted. 

4 Conclusions 

In this paper, we have discussed DIS heavy-flavour production at HERA at NLO. The effect on 
the charm production cross section due to the uncertainties in the factorization/renormalization 
scale n and the charm quark mass has been studied as a function of x and Q^. At small 
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Figure 15: Expected ratio F2^/ F2'^ for a luminosity of 500 pb~^ . The points show the predictions 
according to the calculations of section 2.1 and 2.6. The inner (outer) error represents the 
statistical and the total experimental error. 



X and for < 300 GeV^ the predictions are found to be insensitive to while in this region 

the effect of nic is found to be large. Unfortunately the sensitivity of F2'^{x, Q'^) to the gluon 
density is also restricted to this kinematic region. It has been shown that the contribution of 
Fl has sizeable effects on the charm production cross section at large y. The contribution of 
light quarks to the cross section turned out to be of the order of 5%, nearly independent of x 
and Q^. Results on bottom production have been presented. 

The current experimental situation has been summarized. Based upon this knowledge, 
statistical and experimental systematic uncertainties for large luminosity have been estimated. 
Ultimately an accuracy of 10% in the overall normalization of the cross section may be achieved. 
In the kinematic range where the inclusive measurement F2'^{x, Q^) is found to be sensitive to 
the gluon density, the error of the inclusive measurement will start to be systematics dominated 
with a luminosity of 50 pb~^. It has been demonstrated that exclusive measurements of the 
charm mesons would disentangle the influence of rric and the gluon density on the charm 
production cross section. An exploration of the kinematic plane, the extraction of the gluon 
density, as well as the question when the charrn quark may be treated as a parton, will require a 
luminosity of 500 pb~^. Finally the ratio of F2^/F2^{x, Q"^) has been investigated by performing 
an exclusive D meson analysis. With a luminosity of 500 pb~^ the predicted statistics is found 
to be sufficient to make a detailed study of the x and Q'^ dependence of this ratio. 
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